Substrate recognition and selectivity in SARS-CoV-2 main protease: Unveiling the role of subsite interactions through dynamical nonequilibrium molecular dynamics simulations (preprint)

The main protease (Mpro) of the SARS-CoV-2 coronavirus employs a cysteine-histidine dyad in its active site to catalyse hydrolysis of the viral polyproteins. It is well established that binding of the substrate P1-Gln in the S1 subsite of Mpro active site is crucial for catalysis and this interaction has been employed to inform inhibitor design; however, how Mpro dynamically recognises and responds to substrate binding remains difficult to probe by experimental methods. We thus employed the dynamical nonequilibrium molecular dynamics (D-NEMD) approach to probe the response of Mpro to systematic substrate variations. The results emphasise the importance of P1-Gln for initiating a productive enzymatic reaction. Specifically, substituting P1-Gln with alanine disrupts the conformations of the Cys145 and His41 dyad, causing Cys145 to transition from the productive gauche conformation to the non-productive trans conformation. Importantly, our findings indicate that Mpro exhibits dynamic responses to substrate binding and likely to substrate-mimicking inhibitors within each of the S4-S2' subsites. The results inform on the substrate selectivity requirements and shed light on the observed variations in hydrolytic efficiencies of Mpro towards different substrates. Some interactions between substrate residues and enzyme subsites involve more induced fit than others, implying that differences in functional group flexibility may optimise the binding of a substrate or inhibitor in a particular subsite.

Allosteric modulation by the fatty acid site in the glycosylated SARS-CoV-2 spike (preprint)

The trimeric spike protein plays an essential role in the SARS-CoV-2 virus lifecycle, facilitating virus entry through binding to the cellular receptor angiotensin-converting enzyme 2 (ACE2) and mediating viral and host membrane fusion. The SARS-CoV-2 spike contains an allosteric fatty acid (FA) binding site at the interface between two neighbouring receptor-binding domains. This site, also found in some other coronaviruses, binds free fatty acids such as linoleic and oleic acid, and other small molecules. Understanding allostery and how this site modulates the behaviour of different regions in this protein could potentiate the development of promising alternative strategies for new coronavirus therapies. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to investigate allosteric effects and identify the communication pathways in the fully glycosylated spike in the original SARS-CoV-2 ancestral variant. The results reveal the allosteric networks that connect the FA site to important functional regions of the protein, including some more than 40 [A] away. These regions include the receptor binding motif, an antigenic supersite in the N-terminal domain, the furin cleavage site, the regions surrounding the fusion peptide and a second allosteric site known to bind heme and biliverdin. The networks identified here highlight the complexity of the allosteric modulation in this protein and reveal a striking and unexpected connection between different allosteric sites. Notably, 65% of amino acid substitutions, deletions and insertions in the Alpha, Beta, Delta, Gamma and Omicron variants map onto or close to the identified allosteric pathways.

Antibodies generated in vitro and in vivo elucidate design of a thermostable ADDomer COVID-19 nasal nanoparticle vaccine (preprint)

COVID-19 continues to damage populations, communities and economies worldwide. Vaccines have reduced COVID-19-related hospitalisations and deaths, primarily in developed countries. Persisting infection rates, and highly transmissible SARS-CoV-2 Variants of Concern (VOCs) causing repeat and breakthrough infections, underscore the ongoing need for new treatments to achieve a global solution. Based on ADDomer, a self-assembling protein nanoparticle scaffold, we created ADDoCoV, a thermostable COVID-19 candidate vaccine displaying multiple copies of a SARS-CoV-2 receptor binding motif (RBM)-derived epitope. In vitro generated neutralising nanobodies combined with molecular dynamics (MD) simulations and electron cryo-microscopy (cryo-EM) established authenticity and accessibility of the epitopes displayed. A Gigabody comprising multimerized nanobodies prevented SARS-CoV-2 virion attachment with picomolar EC50. Antibodies generated by immunising mice cross-reacted with VOCs including Delta and Omicron. Our study elucidates nasal administration of ADDomer-based nanoparticles for active and passive immunisation against SARS-CoV-2 and provides a blueprint for designing nanoparticle reagents to combat respiratory viral infections.

Dynamical nonequilibrium molecular dynamics simulations identify allosteric sites and positions associated with drug resistance in the SARS-CoV-2 main protease (preprint)

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalysing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behaviour of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

Molecular dynamics of spike variants in the locked conformation: RBD interfaces, fatty acid binding and furin cleavage sites. (preprint)

Since December 2019 the SARS-CoV-2 virus has infected billions of people around the world and caused millions of deaths. The ability for this RNA virus to mutate has produced variants that have been responsible for waves of infections across the globe. The spike protein on the surface of the SARS-CoV-2 virion is responsible for cell entry in the infection process. Here we have studied the spike proteins from the Original, Alpha (B.1.1.7), Delta (B1.617.2), Delta-plus (B1.617.2-AY1), Omicron BA.1 and Omicron BA.2 variants. Using models built from cryo-EM structures with linoleate bound (6BZ5.pdb) and the N-terminal domain from 7JJI.pdb, each is built from the first residue, with missing loops modelled and 45 disulphides per trimer. Each spike variant was modified from the same Original model framework to maximise comparability. Three replicate, 200 ns atomistic molecular dynamics simulations were performed for each case. (These data also provide the basis for further, non-equilibrium molecular dynamics simulations, published elsewhere.) The analysis of our equilibrium molecular dynamics reveals that sequence variation at the closed receptor binding domain interface particularly for Omicron BA.2 has implications for the avidity of the locked conformation, with potential effects on Omicron BA.1 and Delta-plus. Linoleate binding has a mildly stabilizing effect on furin cleavage site motions in the Original and Alpha variants, but has no effect in Delta, Delta-plus and slightly increases motions at this site for Omicron BA.1, but not BA.2, under these simulation conditions.

SARS-CoV-2 spike variants differ in their allosteric response to linoleic acid (preprint)

The SARS-CoV-2 spike protein contains a fatty acid binding site, also found in some other coronaviruses (e.g. SARS-CoV), which binds linoleic acid and is functionally important. When occupied by linoleic acid, it reduces infectivity, by "locking" the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the response of spike variants to linoleic acid removal. These simulations show that the fatty acid site is coupled to functional regions of the protein, some of them far from the site (e.g. in the receptor-binding motif, N-terminal domain, the furin cleavage site located in position 679-685 and the fusion peptide-surrounding regions) and identify the allosteric networks involved in these connections. Comparison of the response of the original ("Wuhan") spike with four variants: Alpha, Delta, Delta plus and Omicron BA.1 show that the variants differ significantly in their response to linoleic acid removal. The allosteric connections to the fatty acid site on Alpha are generally similar to the original protein, except for the receptor-binding motif and S71-R78 region which show a weaker link to the FA site. In contrast, Omicron is the most affected variant exhibiting significant differences in the receptor-binding motif, N-terminal domain, V622-L629 and the furin cleavage site. These differences in allosteric modulation may be of functional relevance, e.g. in differences in transmissibility and virulence. Experimental comparison of the effects of linoleic acid on different variants is warranted.

A Functional Interaction Between Y674-R685 Region of the SARS-CoV-2 Spike Protein and the Human α7 Nicotinic Receptor (preprint)

The α7 nicotinic acetylcholine receptor (nAChR) is present in neuronal and non-neuronal cells and has anti-inflammatory actions. Molecular dynamics simulations suggested that α7 nAChR interacts with a region of the SARS-CoV-2 spike protein (S), and a potential contribution of nAChRs to COVID-19 pathophysiology has been proposed. We applied whole-cell and single-channel recordings to determine whether the Y674-R685 region of the S protein can directly affect α7 nAChR function. The S fragment exerts a dual effect, acting as a low-efficacy agonist and a non-competitive inhibitor. It activates α7 nAChRs, in line with our previous molecular dynamics simulations showing favorable binding of this accessible region of the S protein to the nAChR agonist binding site. However, activation requires the presence of positive allosteric modulators that enhance open probability, indicating very low efficacy. The main effect of the S fragment on α7 nAChR is a negative modulation, which is evidenced by a profound decrease in the durations of channel openings and activation episodes and in the amplitude of macroscopic responses elicited by ACh. Our study identifies a novel functional interaction between α7 nAChR and a region of the S protein, thus providing molecular foundations for exploring the involvement of nAChRs in COVID-19 pathophysiology.

A Functional Interaction Between the SARS-CoV-2 Spike Protein and the Human α7 Nicotinic Receptor (preprint)

Coronavirus disease 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Infection relies on the binding of the viral spike protein (S) to angiotensin-converting enzyme 2 in host cells. Regions of the S protein have been suggested to interact with nicotinic acetylcholine receptors (nAChRs), and a potential contribution of nAChRs to COVID-19 pathophysiology has been proposed. α7 nAChR is an interesting candidate target since it is present in neuronal and non-neuronal cells, including immune cells, and has anti-inflammatory actions. We here identified a novel direct functional interaction between the α7 nAChR and a region of the S protein (Y674-R685). The S fragment exerts a dual effect, acting as a low-efficacy agonist and a non-competitive inhibitor. It activates the α7 nAChR, in line with our previous molecular dynamics simulations showing favorable binding of this accessible region of the S protein to the nAChR agonist binding pocket. However, activation requires the presence of positive allosteric modulators that enhance channel opening probability, indicating very low activation efficacy. The S fragment exerts an additional inhibition, which may be the predominant effect on α7 responses. This study reveals a functional interaction between α7 and the Y674-R685 fragment of the S protein, which opens doors for exploring the involvement of nAChRs in COVID-19 pathophysiology.

Discovery of SARS-CoV-2 Mpro Peptide Inhibitors from Modelling Substrate and Ligand Binding (preprint)

The main protease (Mpro) of SARS-CoV-2 is central to its viral lifecycle and is a promising drug target, but little is known concerning structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of classical molecular mechanics and quantum mechanical techniques, including automated docking, molecular dynamics (MD) simulations, linear-scaling DFT, QM/MM, and interactive MD in virtual reality, to investigate the molecular features underlying recognition of the natural Mpro substrates. Analyses of the subsite interactions of modelled 11-residue cleavage site peptides, ligands from high-throughput crystallography, and designed covalently binding inhibitors were performed. Modelling studies reveal remarkable conservation of hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular at the P2/S2 sites. The binding modes of the natural substrates, together with extensive interaction analyses of inhibitor and fragment binding to Mpro, reveal new opportunities for inhibition. Building on our initial Mpro-substrate models, computational mutagenesis scanning was employed to design peptides with improved affinity and which inhibit Mpro competitively. The combined results provide new insight useful for the development of Mpro inhibitors.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour (preprint)

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. The results also show significant differences in behaviour between clinical variants of the spike: e.g. the D614G mutation shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Structure and mechanism of SARS-CoV-2 Spike N679-V687 deletion variant elucidate cell-type specific evolution of viral fitness (preprint)

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants which is of particular concern due to their potential for increased transmissibility and pathology. In addition to this entrenched variant diversity in circulation, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriS{Delta} variant, originally identified as a viral subpopulation by passaging SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike glycoprotein an eight amino-acid deletion encompassing the furin recognition motif and S1/S2 cleavage site. Here, we analyzed the structure, function and molecular dynamics of this variant spike, providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity, allowing the virus to probe diverse trajectories in distinct cell types to evolve viral fitness. TeaserSARS-CoV-2 can exploit different cell types to diversify and evolve virus variants distinct in infectivity and structure.

Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors and suggest subtype specificity (preprint)

Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this nicotinic hypothesis. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human 4{beta}2 and 7 subtypes and the muscle-like {beta}{gamma}d receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptors binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the 4{beta}2 and 7 complexes and more compact ones when bound to the muscle-like receptor. In the 4{beta}2 and {beta}{gamma}d complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the 7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.